Depletion forces between non-spherical objects 
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We extend the insertion approach for calculating depletion potentials to the case of non-spherical 
solutes. Instead of a brute-force calculation we suggest to employ the recently developed curvature 
expansion of density profiles close to complexly shaped walls. The approximations introduced in 
the calculation by the use of the curvature expansion and of weight functions for non-spherical 
objects can be tested independently. As an application for our approach we calculate and discuss 
the depletion potential between two hard oblate ellipsoids in a solvent of hard spheres. For this 
system we calculate the entropic force and torque acting on the objects. 
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I. INTRODUCTION 

If macromolecules such as colloids are immersed in a 
solvent of smaller particles, the difference in size makes it 
useful to describe this mixture in terms of effective inter- 
actions by integrating out the degrees of freedom of the 
solvent. The resulting interactions between particles of 
the remaining larger component are often referred to as 
depletion interactions. Such depletion forces have been 
studied in detail both theoretically Q, H 0, 0, IE E| and 
experimentally 0, H ■ Most theoretical approaches are 
based on a brute-force approach in which the solutes are 
frozen in a given configuration and the reby turned into an 
external field for the solvent 0,0 El El El From the 
inhomogeneous structure of the solvent in the external 
field due to two fixed solute particles one can calculate 
the solvent mediated effective force acting on the solutes. 
This brute-force approach turns out to be very time con- 
suming because for each separation and orientation, for 
which one wants to determine the depletion force, the 
inhomogeneous solvent distribution has to be calculated 
anew. 

The insertion approach to depletion potentials dif- 
fers in character in that there only one solute particle is 
fixed. The advantage is that the calculation of the sol- 
vent density profile in the external field of a single solute 
is usually much simpler and hence computational less 
demanding. The potential of the depletion force can be 
determined from the solvent density profile close to one 
solute by inserting the second solute into the system us- 
ing the potential distribution theorem |13| . However, for 
the insertion step a theoretical description of a mixture 
consisting of solute and solvent particles is required. 

So far, most studies of depletion forces have been fo- 
cused on rather simple geometries, such as the force be- 
tween a big sphere and a planar wall or between two big 
spheres in a solvent of small spheres. For these sym- 
metric systems the depletion potential depends on the 
sphere-wall or the sphere-sphere separation as the only 
parameter characterizing the configuration. In addition, 
for these geometries the use of the aforementioned inser- 
tion approach is facilitated by the availability of reliable 
density functional theories for hard-sphere mixtures. 



In colloidal mixtures one generally encounters more 
complex particle shapes both of the solute and solvent 
particles. However, the understanding of depletion po- 
tentials for non-spherical objects is still rudimentary. 
The corresponding calculations are much more challeng- 
ing because in these cases the depletion potential depends 
not only on the separation between the solutes but also 
on their relative orientation. 

A first extension beyond the mixture of spheres is the 
case of spherical solutes immersed in a solvent of non- 
spherical particles. Depletion agents such as thin rods 
[ll E El E3 or infinitely thin platelets El HI can 
generate big depletion effects even at rather low solvent 
concentrations • The depletion potential in these cases 
can be calculated in the limit of small solvent densities so 
that correlations among the depletion agents are small. 
Similar in spirit is also the calculation of the depletion 
force between s phe res in a solvent of a liquid-crystal in its 
nematic phase [2l|, [22 , for which the strong correlations 
between particles of the liquid crystal are taken into ac- 
count effectively by reducing their orientational degrees 
of freedom. The strength of the depletion interaction is 
then estimated by excluded volume calculations, follow- 
ing the ideas of Asakura and Oosawa 0, • 

A second, more complicated situation is the one stud- 
ied in Ref. |23|. where one spherocylinder immersed in a 
spherical solvent close to a planar wall was considered. 
These calculations employed the insertion approach. The 
resulting depletion potential depends not only on the sep- 
aration of the solute from the wall but also on its orienta- 
tion. Hence in addition to the depletion force an entropic 
torque acts on the solute. 

There are recent studies of depletion forces between 
two spherocylinders in a solvent of spheres which over- 
come the Asakura- Oosawa approximation 0. Q Using 
a three-dimensional integral-equation theory, Kinoshita 
[Til Il2| showed that the corresponding depletion poten- 
tial displays a rich behavior and depends sensitively on 
the path along which the spherocylinders approach each 
other. Similar findings were reported in a simulation 
study, in which the depletion potential was determined 
by the acceptance ratio method |24j . 

Here we extend the DFT insertion approach to deple- 
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tion potentials [6| to the case of non-spherical objects. In 
Sec. [HI we recall the basic theory and highlight how the 
geometry of the solutes can be accounted for. We test the 
new elements of the theory in Sec. IIIII As an application 
we discuss the depletion potential between two ellipsoids 
in Sec. IIVI We conclude in Sec. El 

II. THEORY 

We follow the versatile and successful approach to cal- 
culating depletion potentials between two objects, a and 
b, immersed in a solvent within the framework of den- 
sity function theory (DFT) [||, which is referred to as 
the insertion approach. To this end we fix one of the ob- 
jects, say a, at the origin at a given orientation so that 
it acts as an external potential for the solvent particles. 
In response to this external potential, the solvent parti- 
cles aquire an inhomogeneous equilibrium number den- 
sity distribution p s (v). In the fluid phase, the case we 
are interested in here, p s (r) shares the spatial symmetry 
with that of object a. If object a is a sphere, the den- 
sity distribution p s (r) possesses also spherical symmetry. 
In a second step we insert the second object, denoted 
by 6, into the inhomogeneous solvent at position r and 
with relative orientation oj. As a result of this insertion 
the grand potential Q(r, u>) of the system changes. The 
depletion potential is given by 

W(t,u) = fi(r,w) -fi(r oo,u;), (1) 

which can be re-written in terms of the one-body direct 
correlation function cj; (r, w) = —pSJ^ex/Spbir^u}) 0: 

0W(t,u) = 4 1} (r -f oo, w) - 4 1} ('.«)■ (2) 

For the numerical calculation of c[ (r, uj) two challenges 
have to be overcome: (i) the accurate calculation of the 
density profile p s (r) of solvent particles around an object 
of complex shape, and (ii) the insertion of a non-spherical 
object into an inhomogeneous solvent of spheres. Since 
both of these problems require non-standard approaches, 
in the following we shall pay special attention to them. 

Although our approach is flexible and can treat also 
soft solvent-solvent and solute-solvent interactions, in the 
following we shall restrict our considerations to the case 
of hard-core interactions. The solvent is represented by a 
hard-sphere fluid characterized by its radius R and bulk 
density p s or bulk packing fraction rj s — 4-7ri? 3 p s /3. The 
solute-solvent interaction is infinitely repulsive in the case 
of overlap and zero otherwise. 

A. Density profiles and weighted densities 

1. Curvature expansion 

Concerning the first issue of calculating the density 
profile around non-spherical objects we apply the re- 



cently suggested and successfully tested curvature expan- 
sion of the density profile |2j| ■ To this end we introduce 
normal coordinates. Any point r outside the fixed object 
a can be reached from R as the point closest to r on the 
parallel surface of the object, where the density profile 
p a (r) vanishes discontinuously. This particular surface is 
special to the case of hard-core solute-solvent interaction; 
however, it is possible to employ any parallel surface as 
long as all the calculations are in line with this definition 
of R. 

For the vector connecting the points R and r one has 
R r = un(R), where n(R) is the unit vector normal 
to the parallel surface at point R and u is the normal 
distance. At point R the parallel surface exhibits two 
principle radii of curvature R\ and i?2 leading to the 
dimensionless mean and Gaussian curvatures H(R) = 
(R/Ri + R/R 2 )/2 and K(R) = R 2 /(R 1 R 2 ), respectively. 

As an ansatz for the density profile we employ |25| 

Ps (r) = pf («) + H(K)pf («) + K(R)p?(u) 
+H(R) 2 pf(u)+H(-Rfp? 3 (u) 

+H(R)K(K)p™(u)..., (3) 

which factorizes the local geometry of object a at posi- 
tion R, specified by the local mean and Gaussian curva- 
ture H(R) and -ftT(R), respectively, and the structure 
of the solvent via the coefficient functions pj(it) with 
£ = P, H, K, . . .. This separation of geometry and the 
coefficient functions allows one to infer pf (it) from sim- 
ple geometries with high symmetry, such as a fluid close 
to planar, spherical, or cylindrical walls. The functions 
pl{u) have been determined in Ref. [25j. There the direct 
comparison between a density profile predicted by the 
curvature expansion, based on these coefficient functions 
and obtained within DFT, with results from Monte-Carlo 
simulations has demonstrated the high accuracy of this 
approach. Note that the curvature expansion of the den- 
sity profile, Eq. @, implicitly assumes that the curvature 
of the fixed object varies smoothly on its the surface. A 
sharp edge would be problematic because the mean cur- 
vature jumps from zero to a non-vanishing value at the 
edge. 

2. Contact density on spherocylinders 

On the surface of a spherocylinder both the mean and 
the Gaussian curvature vary discontinuously where the 
spherical cap meets the cylinder. Also this discontinuity 
of the curvatures cannot be captured fully by the curva- 
ture expansion. 

In order to analyze the reliability of the curvature ex- 
pansion we have performed a Monte Carlo simulation of 
a hard-sphere fluid exposed to a spherocylinder. As pa- 
rameters we have chosen rj s = 0.3 for the packing fraction 
of the fluid, L — 1QR for the length of the cylinder, and 
R sc = 4i? for the radius of the parallel surface at which 
the density profile discontinuously drops to zero. In this 
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o MC simulations 
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x=(L+R sc ji)/2 




FIG. 1: The contact density p c (x) of a fluid of hard spheres 
with radius R and packing fraction n s = 0.3 at a hard sphero- 
cylinder with length L — 10R and radius R c = 3R so that the 
contact density occurs at R sc = 4R. The path on the parallel 
surface of contact of the spherocylinder is parameterized by 
x, as is indicated by the dotted line in the inset. Note that 
the inset is not drawn to scale. According to the curvature ex- 
pansion the contact value p c (x) jumps at x = L/2, where the 
cylindrical and the spherical parts meet, as shown by the full 
line. The symbols denote data from Monte Carlo simulations 
for rj s = 0.3. 



test we focus on the contact density p c {x) as the most 
sensitive quantity, where x parameterizes a path along 
the surface as depicted in the inset of Fig. The den- 
sity profile away from contact decays towards the bulk 
density and we verified that that the effects of the discon- 
tinuity of the curvature decreases with increasing normal 
distance from the spherocylinder. 

The dimensionless curvatures on the cylindrical part 
of the surface are H = R/(2R SC ) and K = 0, while they 
are H = R/R sc and K = (R/R sc ) 2 on the spherical 
caps. Accordingly, the curvature expansion predicts a 
jump in the contact value of the density profile where 
the spherical caps meet the cylinder. The result of the 
curvature expansion is plotted as the full line in Fig. ^ 
Note that in Eq. © only the three terms pf(u), pf(u), 
and pf (u) have a non- vanishing contact value and hence 
contribute to p c (x) pll2f|. 

In contrast to the jump of the contact density between 
constant values as predicted by the curvature expansion, 
in the computer simulations we find a smooth, slightly 
oscillatory transition between the contact densities at 
the cylindrical and the spherical part of the spherocylin- 
der. In Fig. ^ the simulation data for p c {x) are shown 
as symbols. Interestingly, the spatial region of deviation 
between the simulation data and the prediction of the 
curvature expansion is relatively narrow. 



3. Curvature expansion of the free energy density and of 
weighted densities 

Although the density profile p s (r) , as given by Eq. ||3J| , 
is in principle sufficient for the calculations we intend to 
perform, it is numerically more efficient to exploit the 
particular form of fundamental-measure theory (FMT) 
functionals |27| . Within FMT the one-body direct corre- 
lation function is given by 



®w b a , 



(4) 



where $ is the excess free energy density and ^ a = 
d&/dn a depends on p s (r) in a complicated, nonlinear 
way. The convolution product is denoted as <g>. Instead of 
calculating ^ a from p s (r) directly, we argue that ^a(r) 
can be equivalently expanded in terms of powers of H 
and K and therefore can be written as 



(5) 



One can adopt the point of view that curvature expan- 
sions such as those given in Eqs. J2J and JSJ are merely 
approximations of the functions /0 s (r) or \&(r) that take 
the shape of the external potential into account in an effi- 
cient way. Within this line of arguments there is nothing 
special about the density distribution or any other func- 
tion entering the DFT. Even auxiliary functions such as 
the weighted densities n a (r) and ^ a (r) have a curvature 
expansion ensuring that the output of the DFT calcula- 
tion has the form of Eq. JSJ . 

A more systematic, albeit more involved, point of view, 
which we present here only as a sketch, starts with an ap- 
proach similar in spirit to the one that leads to Eq. © 
[25| as a suitable form for the density profile. Here, how- 
ever, we analyze the weighted densities, defined as 



n a (r) = J dr'p s (r')w a (r,r') 



(0) 



close to planar, spherical, and cylindrical walls. Our 
results suggest that analogous to the density profile 
[Eq. ©] also the weighted densities can be expanded 
in terms of powers of the curvatures H and K . As in 
Ref. |25j we can determine uniquely the coefficient func- 
tions ni for £ = P,H, K, H 2 , HK, H 3 up to third order 
in the inverse of the radii of curvatures inferred from 
planar, spherical, and cylindrical geometries. In order 
to obtain higher order contributions one would have to 
consider more complex wall shapes. The numerical ac- 
curacy of weighted densities calculated at more complex 
walls is, however, unsatisfactory and practically prevents 
a reliable decomposition into coefficient functions. Since 
we are interested in calculating the depletion potential 
between two big non-spherical objects immersed in a sol- 
vent of small spheres, the coefficient functions which we 
have determined are sufficient. For our approach to be 
quantitatively reliable, curvatures of the surface of the 
non-spherical objects should always be sufficiently small. 
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We therefore use as an ansatz for the weighted densi- 
ties, following Ref. 25], 



n a (r)=nP(u) + H{IL)nZ(u)+K(R)nK(i 



(7) 



which we can insert into the free-energy density <f> or its 
derivative with respect to n a , i.e., ty a . Both $ and fy a 
are highly nonlinear functions of n a . However, by Taylor 
expanding $ or ty a into powers of n a and by re-arranging 
terms one can see immediately that the curvature expan- 
sions of ^ a follow directly from Eq. JjJ. 



B. Insertion of non-spherical objects 

We now turn to the second part of the present problem, 
i.e., the calculation of the insertion free energy [Eq. J5J] 
of the non-spherical object b. Note that the change in 
the grand potential due to the insertion of object b into a 
homogeneous bulk fluid at r — > oo can be described by us- 
ing the morphometric approach [2^, [2(| . In the bulk, the 
change in grand potential of the system cannot depend 
on the orientation of the inserted object, which simplifies 
the problem somewhat. Furthermore, it was shown that 
the problems even simplify further due to the separation 
of the geometry and the shape independent thermody- 
namical coefficients. In the morphometric approach, the 
insertion free energy of object & in a bulk fluid can be 
written as El El 
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(r -> oo, u) = pV b + oA b + nC b + RX bl (8) 



where p, a, n, and R are the pressure, the planar wall 
surface tension, and two bending rigidities, respectively, 
which depend on the state of the bulk fluid and the inter- 
action between the fluid and object b. These coefficients 
can be obtained in simple geometries. The corresponding 
geometrical measures V b , A b , C b , and X b describing the 
shape of object b are the volume, the surface area, and 
the integrated (over the surface area) mean and Gaussian 
curvature, respectively. 

In order to calculate the convolutions in Eq. the 
weight functions w b a for a non-spherical object are re- 
quired. For this problem we employ Rosenfeld's formula- 
tion of fundamental measure theory generalized to con- 
vex hard bodies [28l 129) . There are four scalar weight 
functions: 



t4(r)=e(|r-R 6 (M)|), 
which defines the volume V b of object b, 
w b 2 (r)=6(r-R b (9,<t>)), 
which defines the surface area A b of b, 



H(r)w\ 
4tt ' 



(9) 



(10) 



(11) 



which defines the integrated (over the surface) mean cur- 
vature C b of 6, and 



K(r)w b 2 

47T 



(12) 



which defines the integrated (over the surface) Gaussian 
curvature or Euler characteristics X b of object b. Be- 
sides the scalar weight functions, which represent the 
geometrical properties of &, there are two additional 
vector-like weight functions which are required for the de- 
convolution of the Mayer-/ function describing the inter- 
action between non-spherical particles. The vector-like 
weight functions are given by 

w 2 (r) = -Vio|(r) = n b (r)8(v - R b (9, 0)), (13) 

where n b (r) is the unit vector of the surface normal at 
point r, and 

7J(r)w 2 (r) 



wi(r) 



4tt 



(14) 



In a bulk system the vector-like weighted densities n 2 (r) 
and ni(r) vanish. These weight functions have been em- 
ployed successfully in the calculation of the depletion po- 
tential between a hard spherocylinder and a planar hard 
wall (23 - 

With Eqs. J2J and (@J we can now calculate the deple- 
tion potential W(r, u>). 



C. Force and torque 

From the knowledge of the depletion potential W(r, u>) 
it is possible to determine the entropic force and the en- 
tropic torque |23| acting on object b with orientation u> 
at a given position r. If object b is translated by an in- 
finitesimal vector <5r, while keeping its orientation fixed, 
the depletion potential changes by SW = — F • Sr, which 
defines the depletion force 

F(r,w) = -— W(t,u). (15) 

The torque acting on the object b can be calculated by 
a similar consideration, rotating object b by an infinites- 
imal angle Suj. The direction of Su) is parallel to the axis 
of rotation and its modulus specifies the angle of rotation. 
If one keeps the center of b fixed at r and performs a ro- 
tation by Su>, the depletion potential changes by 5W, so 
that SW = — M • Suj. Therefore we can write the entropic 
torque as |23| 

M(t,u)=~W(t,u). (16) 



III. TEST: DEPLETION POTENTIALS 
BETWEEN ELLIPSOIDS AND SPHERES 

Before we apply the above formalism to the calculation 
of the depletion potential between two non-spherical ob- 
jects, we perform a test that enables us to estimate the 
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FIG. 2: Depletion interaction between an oblate (a) [pro- 
late (b)] ellipsoid with half-axes (10, 10, 4)i? [(4,4, 10) R] and 
a sphere with R a — 10R. The center of the inserted sphere 
approaches the north pole of the fixed ellipsoid normal to the 
surface. Both particles touch if u — R a . The symbols denote 
results obtained by approach (1), for which we employ FMT 
for convex objects, and the lines denote corresponding results 
from approach (2), for which we use the curvature expansion. 
As both approaches involve approximations of very different 
nature, the excellent agreement is very likely to be due to 
the fact that the systematic error is very small in both ap- 
proaches. We verified this observation also for ellipsoids with 
different half-axes. For these data the depletion agent is a 
fluid of hard spheres with radius R modeled via the White- 
Bear version of FMT. The bulk packing fraction is r/ a = 0.3. 
Note that the insets here and in the following figures are not 
drawn to scale. 



errors introduced into the numerical calculation through 
the approximations we have made. One source of error 
is the use of the curvature expansion of the density pro- 
file p s (r) [Eq. (PJ] and of the derivatives ^(r) of the 
excess free-energy density [Eq. (J5J]. It introduces an ap- 
proximation because we have to truncate the expansion 
after the third order in the inverse radii of curvature. 
Another source of error is the application of the FMT 
weight functions for non-spherical objects, which intro- 
duces a different approximation. The fact that these two 
approximations are very distinct in nature allows us to 
perform a stringent numerical test. 

To this end we calculate the depletion potential be- 
tween one big sphere and one big ellipsoid in a solvent of 
small spheres. In an exact treatment the depletion po- 
tential between these two objects would depend only on 
their relative position and orientation. In the numerical 
implementation of the insertion method we can follow 
two different routes which make independent use of the 
different approximations. We can choose to fix either the 
sphere and insert the ellipsoid or fix the ellipsoid and in- 
sert the sphere. By the choice of particle fixed at the 
origin we decide about the symmetry of the external po- 



tential and how to calculate the density distribution of 
the solvent. 

If we choose to fix the big sphere the calculation of 
the density profile p s (r) — psix) is straightforward and, 
using the spherical symmetry of the problem, it was es- 
tablished that the results agree extremely well with, e.g., 
Monte Carlo simulations. Along this route, the main ap- 
proximation for the calculation of W(r, ui) stems from 
using the weight functions for the ellipsoid. 

Note that the resulting depletion potential W(r, u>) de- 
pends on both the orientation u) of the ellipsoid relative 
to the vector r connecting the centers of the ellipsoid and 
the sphere and the distance r between the ellipsoid and 
the sphere. For this test, however, we fix the orienta- 
tion and consider the approach between the ellipsoid and 
the sphere along the surface normal of the sphere for the 
chosen orientation of the ellipsoid. In Fig. [5] we show 
the depletion potential between a big sphere, denoted as 
object a, with radius R a = 10 R and an oblate ellipsoid 
with half-axes (10, 10, A)R (a), and between a sphere with 
radius R a — 1QR and a prolate ellipsoid with half axes 
(4, 4, 1G)R (b) . In both cases the solvent is a fluid of small 
spheres with a packing fraction rj s = 0.3, which we model 
by the White Bear version of FMT [3(l|3l|]. The insets in 
Fig. depict the orientation between the ellipsoids and 
the sphere chosen in the calculation. The symbols denote 
the results obtained via the first route, corresponding to 
a fixed sphere. 

Along the second route, we fix the ellipsoid at the ori- 
gin so that it acts as an external potential for the sol- 
vent spheres. Now we employ the curvature expansion 
[Eq. in order to evaluate the derivatives of the excess 
free energy density ^ a (r). The weight functions we need 
in order to describe the insertion of the big sphere in 
Eq. J3J are well tested and are known to be accurate |g] . 
The results for the depletion potentials along the same 
paths between the sphere and the oblate and prolate el- 
lipsoid obtained from this route are shown in Fig. [5] as 
lines. 

We find that the results obtained from both routes 
agree extremely well, which provides confidence in the 
reliability of the approximations and the numerical ap- 
proach. Only in the case of the prolate ellipsoid, for which 
the curvatures are considerably higher than in the case of 
the oblate ellipsoid, we find some deviations between the 
two routes close to the first potential barrier. However, 
these deviations are very small. 

Close to contact between the sphere and the ellipsoids, 
one can employ arguments based on considerations about 
the overlap of excluded volumes [jj, I4 l| ■ In the case of 
high curvature one expects that the contact value of the 
depletion potential is considerably reduced compared to 
cases of low curvature. This expectation is confirmed by 
our results. In addition to the contact value, we find 
that the amplitude of the oscillations of the depletion 
potential for the oblate ellipsoid [Fig.^a)] is larger than 
the one for the prolate ellipsoid [Fig. Hfo)]. 

We conclude from the results of this test that both the 
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FIG. 3: Depletion interaction between two oblate ellipsoids 
with half axes (10, 10, 4)R. The inserted particle approaches 
the north pole of the first ellipsoid perpendicular to its sur- 
face. 7 denotes the angle between the z-axis and the large 
half-axis of the inserted ellipsoid (see the inset). If both el- 
lipsoids are aligned (as drawn in the inset), i.e., 7 = 90°, 
the minimal separation of the centers is z m i n — 8R, while 
for 7 = it is Zmin = 14R; z m i n (~/ — 60°) w 9.53i? and 
Zmin (7 = 30°) w 12.57-R. The positions of the extrema and 
the zeros are basically independent of 7. The solvent is a 
fluid of hard spheres with radius R and bulk packing fraction 
r/ s — 0.3 and is modeled via the White-Bear density func- 
tional. 

curvature expansion of the functions ^ <*(*") and the in- 
sertion of a non-spherical particle into an inhomogeneous 
solvent of small spheres work reliably. 



IV. APPLICATION: DEPLETION POTENTIALS 
BETWEEN TWO ELLIPSOIDS 

We can now turn to the calculation of the depletion 
potential between two non-spherical, convex objects. In 
order to accomplish this we have to combine both steps 
mentioned and tested above. First we fix one of the two 
non-spherical objects and thus turn it into an external 
potential for the solvent of small spheres. The structure 
of the resulting inhomogeneous solvent density distribu- 
tion at a given bulk density is captured by the curvature 
expansion given by Eq. J5J . In the second step we employ 
the insertion of a non-spherical object. 

In order to illustrate our approach we calculate the de- 
pletion potential between two hard oblate ellipsoids with 
half-axes (10, 10,4)i? in a solvent of hard spheres with 
radius R and bulk packing fraction rj s = 0.3. The re- 
sulting depletion potential depends on both the relative 
orientation uj of the two ellipsoids and the vector r con- 
necting their centers and thus depends on 6 variables. 



Out of this high-dimensional parameter space we select a 
few examples of paths along which we study the behavior 
of the depletion potential. 

In the first example, the center of the inserted ellipsoid 
approaches the north pole of the fixed particle along the 
surface normal, as shown in the inset of Fig. [3J For this 
path we vary also the relative orientation a; between the 
particles. Due to the symmetry of this configuration the 
orientation between the oblate ellipsoids can be expressed 
in terms of a single angle, which we denote by 7. If 
we choose 7 = 90°, the ellipsoids are parallel and the 
minimal separation of their centers is z m i n = 8R. At 
contact the overlap of excluded volume is larger than 
for any other value of 7. Hence the contact value of 
the depletion potential is most negative. For our choice 
of parameters we find f3W(z m i n ,j — 90°) ~ —12 (see 
Fig. |3J). In addition to the strong attraction close to 
contact, the depletion potential displays a pronounced 
oscillatory structure away from contact. 

This oscillatory structure of the depletion potential re- 
flects mainly the properties of the solvent. The structure 
of the hard-sphere solvent normal to the surface displays 
oscillatory, exponentially decaying packing effects. We 
have shown previously |6j] that beyond the first maximum 
the depletion potential between two spheres or between a 
sphere and a planar wall decays in a closely related, expo- 
nentially damped oscillatory fashion. The decay length 
and the wavelength of the oscillations in the depletion 
potential are exactly the same as those in the decay of 
the bulk pair correlation function. Only the amplitude 
of this decay and the phase of the oscillations depend on 
the shape and the orientation of the two solutes. 

As we decrease the value of 7 from 90° to 0, thereby 
changing the relative orientation from parallel to normal, 
we find that the contact value as well as the potential bar- 
rier at z — z m i n > R and the amplitude of the oscillations 
decrease monotonically; the positions of the extrema and 
of the zeros basically do not vary with 7. 

In the second example the ellipsoids always approach 
each others north poles in a parallel configuration. The 
corresponding path of the centers is a straight line form- 
ing an angle a with the fixed direction of the small axes 
(see the inset of Fig. 0J. In Fig. ^ by construction the 
path for a = is identical to the path with 7 = 90° in 
Fig. El Upon increasing a, at first the resulting depletion 
potential changes only slightly but for angles a > 80° it 
is possible to observe a clear decrease in the wavelength 
of oscillation, which is most pronounced in the case of 
a = 90°. This, however, is not a contradiction to the 
aforementioned universality of the oscillatory decay of 
the depletion potential, because this path cuts through 
the three-dimensional oscillatory structure of the solvent, 
which is organized normal to the surface of the fixed so- 
lute, at the angle a. 

As pointed out earlier, for non-spherical objects it is 
also possible to keep the distance between their centers 
constant and change the relative orientation. From the 
corresponding change in the depletion potential we can 
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FIG. 4: Depletion potential for the same setup as in Fig. El 
with fixed 7 = 90° . The angle a characterizes the straight line 
along which the center of the inserted ellipsoid is moved away 
from its contact position at the north poles (see the inset). 
The abscissa measures the normal (i.e., minimal) distance u 
of the center of the inserted ellipsoid from the surface of the 
fixed one. For small a the depletion potentials almost coincide 
due to the oblate shape (10, 10, 4)i? of the fixed ellipsoid. 



obtain the entropic torque p3j acting on the inserted 
particle. For the same geometrical setup as in Fig. [3] we 
calculate the torque M for center to center separations 
Az/ R = 10, 12, 14, and 16. Due to the symmetry of this 
setup the entropic torque acting on the inserted ellipsoid, 
relative to its center, is parallel to the rotation a;, i.e., 
M(Az, 7) = Af(Az,7)n 7 , with n 7 = cj/uj and 



M(Az,7) 



dW(Az,j) 
$7 



(17) 



The symmetry of the problem leads to M(Az, 7) = for 
7 = and 90°. A positive value of the torque acts on the 
inserted ellipsoid as to increase the angle 7 (rotating it to- 
wards an orientation parallel to the fixed ellipsoid), while 
a negative value of M leads to a decrease of 7 (rotating 
it towards an orientation normal to the fixed ellipsoid). 
Some typical examples of the torque as a function of 7 
for various values of Az are shown in Fig. |SJ For small 
values of Az the amplitude of the torque is largest; how- 
ever, the hard-core interaction prohibits small values of 
7 due to geometrical constraints. For larger separations 
Az between the ellipsoids, the accessible range of values 
of 7 increases until finally the inserted ellipsoid can ro- 
tate freely without encountering the fixed ellipsoid. Due 
to the symmetry of the problem the torque vanishes for 
7 = 90° as well as for 7 = 0, provided this orientation is 
accessible. 

If for any relative orientation of the ellipsoids the min- 
imal distance between the surfaces equals the diameter 
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FIG. 5: Entropic torque between two oblate ellipsoids with 
half-axes (10, 10, 4)R immersed in a solvent of hard spheres 
with radius R as sketched in the inset. The angle 7 param- 
eterizes the rotation of the inserted ellipsoid and Az denotes 
the distance between the centers. The ellipsoids are aligned 
such that the vector connecting their centers runs through the 
north poles of the fixed (lower) ellipsoid. For this setup the en- 
tropic torque M is given by M = -dW/d'y [1|. If Az > 8R, 
the inserted (upper) ellipsoid can be rotated as indicated in 
the inset and experiences an entropic torque. Positive values 
of M mean that the ellipsoid is pushed towards larger val- 
ues of 7. The curves end if the two ellipsoids overlap due to 
the geometrical constraint. Note that the visible discontinu- 
ities in the first derivative of the curves for Az = 12R and 
Az = 14_R are not artifacts. They occur if both ellipsoids 
come so close that just a single small sphere with radius R 
fits in between their surfaces. For Az = 16i? the distance 
between the ellipsoids is large enough so that this effect does 
not occur. 



2R of a solvent spheres, the torque exhibits a cusp, as can 
be seen in the cases of Az > 10R. For Az = 1QR this ori- 
entation occurs at 7 = 90°. Accordingly, for sufficiently 
large Az no such cusps occur. 

All these feature of the torque (see Fig. [SJ arc similar 
to those of the entropic torque acting on a spherocylinder 
close to a planar wall p3 |. 

From these results, one can speculate how two freely 
floating ellipsoids most likely would approach each other. 
Similar to the case of a spherocylinder close to a planar 
wall |2jj we find that the ellipsoids prefer to be parallel 
once they touch each other. This configuration corre- 
sponds to 7 = 90° in Fig. [21 In order to be able to 
overcome potential barriers while approaching, it is eas- 
ier for the ellipsoids to approach with a relative orien- 
tation different from the parallel configuration. By first 
seeking contact in regions of high curvature, and then by 
adapting the orientation such that regions of low curva- 
ture get in contact with each other, potential barriers, 
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which have to be overcome, are significantly reduced and 
the minimum of the depletion potential can be reached. 
This interpretation agrees with observations reported in 
Refs. |0, 0, 0| where the depletion potential between 
two spherocylinder was studied. 

V. SUMMARY AND CONCLUSIONS 

We have extended our density functional theory ap- 
proach for calculating depletion potentials to the ef- 
fective interactions between generally shaped convex par- 
ticles with surface curvatures which vary smoothly. As 
in the previously studied geometrically simple case of 
two spherical particles, it is most efficient to carry out the 
calculation in two steps. In the first step we describe the 
structure of the solvent close to a fixed object. We have 
shown that a curvature expansion of the density profile, 
or equivalently of the derivatives \& a of the free energy 
density [Eq. (JSJ] provides an efficient tool. By employing 
the curvature expansion we implicitly assume that the 
mean and Gaussian curvatures of the surface of the so- 
lute vary smoothly across its surface. Discontinuities in 
the curvatures, such as those observed at sharp edges or 
close to the spherical caps of a spherocylinder cannot be 
captured fully by the curvature expansion in its present 
form (see Fig. 

In a second step we insert the second particle into the 
inhomogeneous solvent. To this end we employ an ex- 
tension of fundamental measure theory to non-spherical 
particles [2^, Hj| , which was already used in the study 



of a spherocylinder close to a planar wall [30(. Using 
this approach we have studied the depletion potential 
between one ellipsoid and a big sphere immersed in a 
hard-sphere solvent. In this case we could test the ac- 
curacy of our approach by carrying out the calculations 
in two different ways, which employ the various approx- 
imations of our approach independently. We have found 
excellent agreement between the two routes (see Fig. |2J) 
which provides confidence for the scheme used. Further- 
more, we have studied the depletion potential between 
two equally sized oblate ellipsoids. For this case we have 
illustrated the potential of our approach. From the re- 
sulting depletion potential (see Figs. and ^} one can 
calculate the entropic force acting on the centers of the 
solutes as well as the entropic torque (Fig.EJ. This pro- 
vides a picture for the likely pathway of how two freely 
floating ellipsoids approach each other under the action 
of entropic forces 

Besides the application to colloidal mixtures of non- 
spherical objects and spheres, this approach should prove 
useful to studying biological inspired model key and lock 
systems ^jj, f° r which depletion interactions between 
non-spherical objects and geometrically structured sub- 
strates are considered. The non-spherical objects and the 
substrates display a perfect geometrical match, similar to 
biological macromolecules (key) which can form chemi- 
cal bonds with a cavity (lock) only if they are sufficiently 
close together. The issue is how the key can be guided 
into the lock using a robust, chemically unspecific mech- 
anism. Our present analysis provides first steps towards 
addressing this issue. 
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